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FLUTTER AT MACH 3 OF THERMALLY STRESSED PANELS 
AND COMPARISON WITH THEORY FOR PANELS 


WITH EDGE ROTATIONAL RESTRAINT 

By John L. Shideler, Sidney C. Dixon, 
and Charles P. Shore 
Langley Research Center 

SUMMARY 

An investigation was conducted at a Mach number of 3 in the Langley 9- by 6-foot 
thermal structures tunnel to study effects of thermal stress and buckling on the flutter 
characteristics of initially flat isotropic panels with length-width ratios of 2.5 and 2.9. 

The experimental results are presented in tabular form and are summarized graphically 
in terms of nondimensional flutter parameters. Comparison of experimental results with 
theory that accounts for edge rotational restraint showed good agreement at zero midplane 
stress and over a major portion of the flat panel flutter boundary. The theory was based 
on two-dimensional static aerodynamics. In the region of the transition from the flat 
panel boundary to the postbuckled boundary the agreement becomes poor. However, it is 
believed that accounting for panel initial imperfections in the theory may improve agree- 
ment between experiment and theoretical flutter boundaries in the transition region. 

INTRODUCTION 

The establishment of an adequate design criteria for the prevention of panel flutter 
has proved difficult because of the large number of parameters that affect the flutter 
characteristics of panels. One parameter that appears to be particularly significant is 
midplane compressive stress. The flutter of stressed panels has been the subject of a 
number of theoretical and experimental investigations. (See, for example, refs. 1 to 15.) 
These studies have indicated that considerable discrepancy may exist between theory and 
experiment. Theory predicts that the minimum dynamic pressure for flutter may occur 
at a value of midplane load less than that required for buckling. (See ref. 13.) However, 
existing experimental data show that the minimum value of dynamic pressure occurs at 
the transition from the boundary for flat unbuckled panels to the boundary for buckled 
panels. Most data on stressed panels have been obtained under aerodynamic heating con- 
ditions where thermal stresses have not been accurately known. Thus, the amount of dis- 
agreement between theory and experiment is uncertain. 



In the present investigation, the effect of thermally induced stress on the flutter 
characteristics of flat rectangular panels is again considered with emphasis on those 
areas where discrepancies between theory and experiment exist. Tests were conducted 
at a Mach number of 3 in the Langley 9- by 6-foot thermal structures tunnel to obtain 
experimental data for panels with length- width ratios of 2.5 and 2.9. The experimental 
results presented in tabular form are summarized in terms of nondimensional flutter 
parameters and are compared with theory that accounts for panel edge rotational restraint. 
The theory is a closed-form solution which incorporates two-dimensional static aero- 
dynamics. The average degree of rotational restraint was determined from the measured 
panel natural vibration frequencies. 


SYMBOLS 


The units used for the physical quantities defined in this paper are given both in the 
U.S. Customary Units and in the International System of Units, SI (ref. 16). The appendix 
presents factors relating these two systems of units. 
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panel length (longitudinal direction parallel to airflow) 
width of frame and panel (see fig. 1) 

panel width (lateral directional perpendicular to airflow) 
constant (see eq. (1)) 
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panel flexural stiffness, — -t 

12(l - M 2) 

Young's modulus 
flutter frequency 

natural frequency of nth mode; n = 1,2, 3, 4 

first natural frequency of a simply supported semi-infinite plate, 
panel thickness 
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inplane stress resultant in x-direction (positive in compression), a x h 

inplane stress resultant in y- direction (positive in compression), Oyh 

static differential pressure acting on panel skin (positive pressure deflects 
skin toward airstream) 

free- stream dynamic pressure 

a0 v 

rotational restraint coefficient on leading and trailing edges, 

rotational restraint coefficient on side edges, 

N a 2 

nondimensional loading parameter, x 

tt 2 D 

panel skin temperature 
free- stream stagnation temperature 
average increase of panel skin temperature 
time 

panel coordinate in direction parallel to airflow 
panel coordinate in direction perpendicular to airflow 
coefficient of thermal expansion of panel skin 


0 = v/m 2 - 1 


y mass density of panel skin 

9 X ,6 y rotational spring constants along the leading and trailing edges and along the 

side edges, respectively 

H Poisson's ratio (taken equal to 0.3) 

CT x, CT y midplane stress in x- and y-direction, respectively (compression positive) 
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modified temperature parameter, 
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Subscripts: 


n natural mode shape number (number of half-waves in x- direction) 


T at transition point 

EXPERIMENTAL INVESTIGATION 


Panels 


The panels consisted of flat sheets of 2024-T3 aluminum alloy of various thicknesses 
riveted to thick frames of the same material. The skin and frame were separated by a 
0.031 inch (0.08 cm) strip of fiber glass insulation. In order to reduce initial stresses 
due to mounting, the frames were fastened to a mounting fixture (used to hold the panels 
during tests) before the sheets were riveted. For details of panel construction and 
mounting arrangement, see figures 1 to 4. The panels were 26 inches (66 cm) long and 
10.4 and 9.0 inches (26.4 and 22.8 cm) wide (measured between the center lines of the 
rivets) corresponding to length-width ratios of 2.5 and 2.9, respectively. 
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Figure 1.- Panel construction details (typical of all panels). All dimensions are in inches (cm). 
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L- 64- 2337.1 

Figure 2.- Panel holder in test section as 
viewed from upstream. 

Test Apparatus 

Tunnel .- All tests were 
conducted in the Langley 9- 
by 6-foot thermal structures 
tunnel, a Mach 3 intermittent 
blowdown facility exhausting 
to the atmosphere. A heat 
exchanger is preheated to 
provide stagnation tempera- 
tures up to 660° F (622° K). 
The stagnation pressure can 
be varied from 60 to 200 psia 
(414 to 1379 kN/rn 2 abs). Ad- 
ditional details on the tunnel 
are presented in reference 8. 



(a) Leading and trailing edge detail. 



(b) Side edge detail. 

Figure 4.- Panel mounting arrangement (typical of all panels). 


Panel holder and mounting arrangement .- The panels were mounted in a panel 
holder which extended vertically through the test section (fig. 2). The panel holder has a 
beveled half-wedge leading edge with a recess on the nonbeveled side 29 inches (73.6 cm) 
wide, 30 inches (76 cm) high, and 5 inches (12,7 cm) deep for accommodating test speci- 
mens (fig. 3). Installation of instrumentation in the cavity and instrumentation chamber 
reduces the depth to approximately 3.5 inches (9 cm). Pneumatically operated sliding 
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doors protect test specimens from aerodynamic buffeting and heating during tunnel start 
and shutdown. Aerodynamic fences prevent shock waves emanating from the doors from 
interfering with the airflow over the test specimen. Results of pressure surveys indicate 
that the flow conditions over the exposed surface of a flat panel are essentially free- 
stream conditions (ref. 8). A vent door on the side opposite the recess is used to control 
the pressure inside the cavity behind the test specimen (fig. 3). To improve control of 
the differential pressure, all other openings to the cavity were sealed. 

All panels were mounted flush with the flat surface of the panel holder. The test 
panels and the filler plates were attached with screws to the mounting fixture which had 
been bolted to the panel holder. (See figs. 2 and 4.) 

Instrumentation 

Iron-constantan thermocouples spotwelded to the back of the panel skins at 19 loca- 
tions (see fig. 5) were used to measure panel temperatures. Variable-reluctance-type 
deflectometers were used to detect motion of the panel skin and to measure flutter fre- 
quencies. The deflectometers were located in the cavity approximately 1/4 inch 
(3/5 cm) behind the panel at the three positions indicated in figure 5. In addition, high- 
speed 16- mm motion pictures provided supplementary data on the behavior of the panels. 
Grid lines were painted on the panels for photographic purposes. 
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Figure 5.- Location of panel instrumentation (typical of all panels). All dimensions are in inches (cm). 
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Quick response strain- gage- type pressure transducers were used to measure static 
pressures at various locations on the panel holder and in the cavity behind the panel. 
Stagnation pressures in the test section were obtained from static pressure measurements 
in the tunnel settling chamber. Stagnation temperatures were measured by total tempera- 
ture probes located in the test section. For each test, data from temperature and pres- 
sure transducers were recorded on magnetic tape every twentieth of a second. Deflec- 
tometer readings were recorded on a high-speed oscillograph. 

Test Procedure 

The panels were vibrated at sea-level conditions in the panel holder prior to each 
test with the air- jet shaker described in reference 17. Several panels attached to the 
mounting fixture were also vibrated prior to being installed in the panel holder so that the 
cavity behind the panels was essentially infinite. Comparison of the results indicated lit- 
tle effect on the natural vibration frequencies due to the change in cavity depth. The natu- 
ral frequencies measured in the panel holder prior to each test are presented in table I. 
Also included are the theoretical natural frequencies for clamped and simply supported 
panels calculated from reference 18. Generally, the panel frequencies fell between the 
theoretical results for clamped and simply supported panels, and most frequencies were 
within 10 percent of the value for a clamped panel. No panel was tested if its frequen- 
cies differed from the theoretical value for a clamped panel by more than 25 percent. 

The mode shapes associated with the four frequencies f n recorded during the vibration 
tests consisted of one half-wave in the cross- stream direction and n half-waves in the 
streamwise direction. 

The wind-tunnel tests were conducted at a Mach number of 3, at dynamic pressures 
from 1500 to 5000 psf (72 to 240 kN/m^), and at stagnation temperatures from 300° F to 
500° F (420° K to 530° K). The protective doors on the panel holder were opened after 
the desired test conditions were established and were closed 3 seconds prior to tunnel 
shutdown. The duration of test conditions varied from approximately 10 to 40 seconds. 
The stagnation temperature was essentially constant during each test. The differential 
pressure AP was maintained as near zero as possible by using a manual control. The 
dynamic pressure was constant during the first few seconds of each test, and only thermal 
stresses due to aerodynamic heating varied. However, during the latter part of several 
tests the dynamic pressure was varied in an attempt to obtain as many flutter points as 
possible. The occurrence of flutter was readily determined by monitoring the deflec- 
tometer traces on the high-speed oscillograph during the tests. The usual procedure 
for varying the dynamic pressure was as follows: 

(a) If no flutter had occurred after a predetermined period of time, either the test 
was ended or the dynamic pressure was increased in an attempt to initiate flutter. 
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(b) If flutter had started and stopped, the dynamic pressure was increased in an 
attempt to restart flutter. 

(c) If the panel was still fluttering after a predetermined period of time, the 
dynamic pressure was decreased in an attempt to stop flutter. 

RESULTS AND DISCUSSION 

The results of the investigation are presented in tables II and III in terms of the 
panel and wind-tunnel conditions at the beginning and cessation of flutter. The data tabu- 
lated include the stagnation temperature T^, dynamic pressure q, panel skin tempera- 
ture increase AT, differential pressure AP, and flutter frequency f. Data from tests 
in which no flutter occurred are tabulated at the highest value of the skin temperature 
increase AT attained (these tests were made at a constant dynamic pressure). 

Panel Temperatures 

At the beginning of each test, the panel skin and supporting structure were essen- 
tially at the same temperature. After a panel was exposed to the airstream, the skin 
temperature increased in a manner similar to the typical panel temperature history 
shown in figure 6. The upper 
curve represents the average of 
thermocouples 2 to 9, 12, 13, and 
16 to 19. Individual temperatures 
for these thermocouples were 
within 10° F (5° K) of the average 
value. The lower two curves, 
which represent the average of 
thermocouples 11 and 14, and 1, 

10, and 15, indicate a temperature 
gradient near the panel edges. 

These gradients, which are attrib- 
uted to heat conduction to the sup- 
porting structure, were not 
accounted for in analyzing the test 
data. The average increase in 
temperature of thermocouples 2 to 
9, 12, 13, and 16 to 19 was used as 
the uniform panel skin tempera- 
ture increase AT which is em- 
ployed in subsequent calculations. 
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Figure 6.- Measured panel temperature history for test 12, typical of all tests. 
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Flutter Parameters 


The flutter data obtained in this investigation are summarized in terms of a dimen- 
sionless dynamic pressure flutter parameter and a dimensionless modified temperature 

parameter. Of the quantities in the flutter parameter p which results from 


theoretical analyses employing two-dimensional static aerodynamics such as reference 1, 
only the dynamic pressure q and skin thickness h were varied in these tests. 

Because of the short duration of the tests and the relatively low panel temperatures, 
changes in material properties with temperature were assumed to be negligible. 


Theoretical analyses indicate that the effect of midplane stress on the flutter 


parameter 


s _ N x a 


/ q \1/ q — — 

] - is a function of the loading parameter R x where R x = — = — . 

y3E J h 7T"D 


If N x is assumed to be caused by complete restraint of thermal expansion on all edges 
corresponding to a uniform temperature increase AT, the loading parameter is equiva- 


lent to 12gAT (l * Hi 

7J-2 

this parameter is still useful as a measure of the actual thermal stress. Further, a dif- 
ferential pressure across the panel can also cause midplane stresses. This effect has 
been approximately accounted for in reference 9 by the addition of a pressure term to the 
temperature term to give 


Although thermal expansion was not completely restrained, 


Xfy = 


12(1 + M) 

7r2 



± c 


|AP|/a\ 4 ] 2 ^®\ 
_ E f 


( 1 ) 


This modified temperature parameter \p represents a measure of the midplane loading 

parameter R x . In the pressure term the minus sign applies when a panel is flat because 

a differential pressure causes tension. However, when a panel is buckled, the plus sign 

applies because xj/ is then a measure of buckle depth and both AT and AP tend to 

increase the depth of buckle. The constant c is a function of 3: and is determined 

b 

from the experimental data by the method presented in reference 9. The empirical values 
obtained were 0.60 and 0.28 for the 2.5 and 2.9 length-width ratio panels, respectively. It 
should be emphasized that this differential pressure correction is only an approximation 
and becomes less accurate as the absolute value of AP increases. 


Experimental Flutter Results 

Flutter boundaries.- The flutter boundaries are shown in figure 7 (- = 2.5] and 

fa. \ /a \l/3 a \ b / 

figure 8 - = 2.9 in terms of the flutter parameter |— ] — and the modified tempera- 

\h / \0E/ h 

ture parameter \]s. The open symbols represent flutter- start points (panel flat), whereas 
the open symbols with flags denote that the panel was buckled when flutter started. The 
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O Flutter start, panel flat 
• Flutter stop, panel buckled 

□ Small amplitude flutter start, panel buckled 
■ Small amplitude flutter stop, panel buckled 
3 No flutter, constant dynamic pressure test 
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Figure 7.- Variation of flutter parameter with thermal stress. ^ = 2.5. 
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3 No flutter, dynamic pressure test 
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Figure 8.- Variation of flutter parameter with thermal stress. ^ = 2.9. 
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flutter- stop points (panel buckled) are represented by the solid symbols. The half- solid 
symbols denote the maximum value of i// attained during tests in which no flutter 
occurred. These tests were made at a constant dynamic pressure, as previously noted. 
The curves shown are boundaries which have been faired through the experimental points. 
Each boundary consists of a flat panel portion and a buckled panel portion which intersect 
at a transition point, the minimum point of each boundary. For the panels with length- 

width ratio of 2.5 (fig. 7), the value of ^ ^ the transition point (2.15) is 49 per- 

cent of its value at zero stress (4.35). Beyond the transition point, the boundary slopes 
upward. This behavior is attributed to the gradual increase of the stiffness as the buckle 
depth increases. For the panels with length-width ratio of 2.9 (fig. 8), the value of 

(jtt) ^ h trans *tt° n P° int (2.55) is 54 percent of its zero stress value (4.70). The 

"panel buckled" portion of the boundary slopes upward as did the boundary for ^ = 2.5. 
The general trend of each boundary is similar to previous experimental results. (See, 
for example, refs. 9 and 12.) 

The square symbols at large values of \p in figures 7 and 8 correspond to flutter- 
start and flutter- stop points of small amplitude motion about a buckled shape. This type 
of flutter will be discussed in the next section. 

Flutter behavior .- The flutter observed from high-speed motion pictures appeared 
to be of the traveling wave type. However, in several tests, the motion during the first 
few seconds of flutter appeared to be more characteristic of the standing wave type. Pre- 
dominantly standing wave flutter has been observed for panels with length-width ratios of 
2 or less (refs. 9 and 12), whereas panels with length-width ratios of 4 or more have 
appeared to flutter in a traveling wave (refs. 7, 8, and 10). It follows that panels with 
length- width ratios between 2 and 4 could flutter in either a standing wave or a traveling 
wave or, perhaps, both. The flutter mode shape appeared to have two half-waves in the 
stream direction and one half-wave in the cross- stream direction as did the buckled mode 
shape at all flutter- stop points. This similarity between the flutter mode and buckling 
mode has been observed previously (refs. 7 and 8). 

The flutter amplitude for panels with low stress was found to be at a maximum near 
the trailing edge as is indicated by both theory and experiment (refs. 2 and 10). For 
example, the amplitude and phase relation of the deflectometer records shown in 

figure 9(a) for test 4 ^ = 2.5^ indicate that the maximum amplitude of motion lies near the 

center and between the center and rearward deflectometer locations. However, for panels 
with appreciable thermal stress when flutter started, the largest amplitude was found to 
be in the center as is indicated by the deflectometer traces shown in figure 9(b) for test 5 

^ = 2.5j. As the panel stress approached the buckling stress, the maximum amplitude at 
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flutter start was observed to be 
nearer to the forward deflectom- 
eter than to the rearward deflec- 
tometer as indicated by the deflec- 
tometer amplitude and phase 
relation shown in figure 9(c) for 

test 8 ^ = 2.5^. Theory indicates 

that the maximum amplitude can 
move upstream with an increase in 
stress (ref. 19). This phenomenon 
apparently has not been noted in 
previous experimental work and, 
generally, it was believed that no 
change in mode shape occurred as 
stress increased. (See ref. 13, for 
example.) 
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(a) Maximum amplitude rearward of center (test 4, ~ = 2.5^. 
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Motion pictures of the small 
amplitude flutter represented by 
square symbols in figures 7 and 8 
revealed that the motion appeared 
as a standing wave superimposed 
on the buckle shape, with two half- 
waves perpendicular to the airflow. 
Because of its small amplitude, 
this flutter may not be damaging. 
However, scarcity of data prohibits 
any definite conclusions regarding 
the significance of this flutter. 
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(c) Maximum amplitude forward of center (test 8, | = 2.sj. 

Figure 9.- Sample deflectometer records (plus sign indicates outward 
panel deflection). 


Analytical Results and Comparison With Experiment 

As mentioned previously, the measured natural frequencies generally fell between 
the frequencies calculated for clamped edges and simply supported edges, indicating that 
the panel edge conditions were between those for clamped and simply supported. For 
this reason, experimental flutter results are compared with theoretical flutter results 
calculated from the closed-form solution of reference 1 for simply supported panels and 
with the approximate analysis of reference 2 for clamped panels. It should be noted that 
in reference 2 a cross- stream deflection function was assumed in order to reduce the 
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partial differential equation to an ordinary differential equation. The solution to the 
resulting ordinary differential equation, however, is exact. Clamped-beam modes from 
reference 20 were used as the cross- stream deflection function in the calculations for the 
clamped panel. 

Calculations from these solutions, which use two-dimensional static aerodynamics, 

are compared with experimental data for - = 2.5 and - = 2.9 in figures 10 and 11, 

b b 

respectively. The dynamic pressure parameter is shown as a function of stress for 
simply supported and clamped panels with stress ratios of 1. The values of R x at 
transition R^-p were determined as those values of R^ at which the flutter frequency 
is zero. It should be noted that the curves in figures 10 and 11 are shown as functions of 
a stress ratio and not of stress itself. The value of R x >p is smaller for the simply sup- 
ported panel than for the clamped panel. Thus, at the same value of the abscissa, the 
stress in the panel is much higher if clamped edges are assumed rather than simply sup- 
ported edges. 

The "panel flat* portions of the experimental boundaries (from figs. 7 and 8) are 
included in figures 10 and 11 and are normalized by the value of the modified temperature 
parameter at the transition point For the condition of zero stress, figures 10 and 

11 show that the experimental values of the flutter parameter are bracketed by the 




Figure 10.- Comparison of flat panel portion of faired experi- 
mental flutter boundary from figure 7 with theoretical 
boundary for panels with clamped and simply supported 

edges, f = 2.5; £ = 1.0. 


Figure 11.- Comparison of flat panel portion of faired experi- 
mental flutter boundary from figure 8 with theoretical 
boundary for panels with clamped and simply supported 

edges - 1 = 2 - 9; 5 = i -°- 
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theoretical results for clamped and simply supported edges. Although agreement between 
theory and experiment is good at zero stress, it becomes poor as stresses approach buck- 
ling. For example, the theoretical clamped-panel flutter boundary for - = 2.9 (fig. 11) 

b 

decreases to zero and then increases to the transition value. The theoretical results for 
clamped and simply supported panels presented in figures 10 and 11 suggest a large influ- 
ence due to panel edge rotational restraint. Moreover, the measured natural frequencies, 
which were generally less than theory for panels with clamped edges, indicate that the 
experimental panels were neither clamped nor simply supported. The degree of rota- 
tional restraint which existed in the test panels was estimated for each test by using the 
results of reference 18, which gives panel natural vibration frequencies as a function of 
rotational restraint. The results of reference 18 are given in terms of nondimensional 

a0 v b0 v 

rotational restraint coefficients q x = and q y = where 8 X and 0 y are spring 

constants. Because of the uniformity of edge attachment it was assumed that the panels 
were supported with equal rotational restraint on all edges, that is, 8 X = By. Thus, q y 
is related to q x by q y = q x . The measured first natural frequencies of the test panels 
were used as inputs to calculate corresponding values of rotational restraint coefficient 
by using the numerical results from reference 18. From these values, average rotational 

restraint coefficients were determined to be q x = 85 and q y = 34 ^for -^ = 2 . 5 j and 

q x = 44 and q y = 15 ^for ^ = 2.9^. A few panels, those having the smallest thicknesses, 

had natural frequencies very near the theoretical values for clamped edges. (See, for 
example, tests 1 and 2, table 1(a).) Flutter points for these panels usually fell above the 
faired experimental boundary. Since the boundary was established primarily on the basis 
of tests on panels with finite rotational restraint, the thinner panels that appeared to be 
essentially clamped (indicated by an asterisk in table I) were ignored in establishing the 
average rotational restraint. 

The average of these rotational restraint coefficients and the theory of reference 19, 
which accounts for arbitrary degrees of rotational restraint, were used to calculate theo- 
retical flutter boundaries. The theoretical transition points were again those values of 
R x at which the flutter frequency is zero. 

Experimental boundaries are compared with theoretical flutter boundaries in 
figure 12 ^ = 2.5^ and in figure 13 ^ = 2. 9^ for panels having the calculated finite rota- 
tional restraint. These theoretical boundaries are in reasonable agreement with experi- 
ment not only at zero stress, but also over the major portion of the boundary. Only in 
the region of transition is the agreement poor. Comparison of figures 12 and 13 with 
figures 10 and 11, respectively, indicates that the inclusion of the effects of finite rota- 
tional restraint improves the agreement between theory and experiment. Thus, finite 
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Figure 12.- Comparison of flat panel portion of faired experi- 
mental flutter boundary from figure 7 with theoretical 
boundary for panel with finite edge rotational restraint, 
a N y 

b = 2#5; = i,0; 0 x = 0 y; q X = 85 * 



Figure 13.- Comparison of flat panel portion of faired experi- 
mental flutter boundary from figure 8 with theoretical 
boundary for panel with finite edge rotational restraint. 

5 = 2.9; ^ = 1-0; S x = 6y; Qx=44. 


rotational restraint has an important effect on the flutter behavior of a panel, especially 
in the region of high stress. (Note that the zero dynamic pressure flutter point indicated 
for the panel with clamped edges in figure 11 does not appear in figure 13.) Other factors 
such as length-width ratio and stress ratio also influence the flutter boundary, and more 
accurate definition of these factors may improve the agreement. Inclusion of aerodynamic 
and structural damping in the theoretical calculations can affect the results, but the effects 
of damping become small near the transition point where the greatest need for improve- 
ment exists. 


Another factor, initial imperfections, can have an influence on the flutter behavior 
of stressed panels. For example, a two- mode analysis using two-dimensional quasi- 
steady aerodynamics (ref. 6) has shown that initial imperfections have a large effect on 
the flutter of compressively stressed semi-infinite panels. Thus, initial imperfections 
might also be the cause of the discrepancy between theory and experiment near the tran- 
sition point for the panels of this investigation. Consideration of theoretical and experi- 
mental flutter frequencies also suggests the same cause. The variation of experimental 
flutter- start frequencies with stress is presented in figures 14 and 15 in terms of the 
dimensionless frequency ratio f/f 0 and the dimensionless stress ratio i/ytyrp. Experi- 
mental frequencies in the buckled region are not shown. Theoretical flutter frequencies 
calculated from reference 19 for the two values of a/b are shown for comparison. At 
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Figure 14.- Comparison of theoretical and experimental flutter frequencies. 
F = 2 - 5 = 5 = LO; ®X = e y; <lx = 85. 



low stress the experimental flutter-start 
frequencies are in reasonable agreement 
with theory. At buckling, however, the 
theoretical flutter frequencies go to zero 
as is always true for a linear small 
deflection flutter analysis. The flutter- 
start frequencies did decrease as the 
stress at the flutter- start points 
increased, but these frequencies did not 
approach zero. The pronounced effect 
that initial imperfections have upon the 
natural vibration frequencies of stressed 
beams and plates is illustrated in refer- 
ences 21 and 22. Inclusion of initial 
imperfections in the analysis for 
stressed beams (ref. 21) eliminated the 
zero frequency values at buckling. The 
resulting theoretical natural frequencies 
agreed very well with measured frequen- 
cies. The effects of initial imperfec- 
tions were also considered to contribute 
appreciably to the discrepancies 
between theoretical and experimental 
natural vibration frequencies for a plate 
near the region of buckling. Thus, the 
neglect of initial imperfections in the 
present panel flutter theory may be the 
cause of the discrepancy near transition. 

CONCLUDING REMARKS 

An investigation was conducted at 
a Mach number of 3 in the Langley 9- by 
6 -foot thermal structures tunnel to 
study the effects of thermal stress on 
the flutter characteristics of flat iso- 
tropic panels with length-width ratios of 
2.5 and 2.9. The results revealed that 
the flutter boundaries consisted erf a flat 
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panel portion, a buckled panel portion, and a transition point at their intersection. The 
boundaries are similar to experimental boundaries previously obtained for panels with 
length-width ratios from 0.96 to 10. 

The flat panel flutter boundaries were compared with small deflection theory for 
simply supported and clamped edge conditions. Whereas agreement between theory and 
experiment was good at zero stress, it became poor as stresses approached buckling. 

By incorporating rotational restraint in theory, agreement between experimental and theo- 
retical flutter boundaries was improved greatly; thus, the degree of edge rotational 
restraint has an important effect on panel flutter behavior. The faired experimental 
boundary was in reasonable agreement with theory for elastically restrained panels not 
only at zero stress but also over a major portion of the boundary. In the region of the 
transition point agreement remained poor. 

Comparison of experimental flutter frequencies with theoretical flutter frequencies 
that account for rotational restraint revealed that for the condition of zero stress, the 
experimental frequencies were in good agreement with theory. However, for stresses 
that cause buckling, the experimental frequencies were considerably greater than the zero 
value predicted by the theory. It is believed that the inclusion of panel initial imperfec- 
tions in the theory may improve the prediction of experimental and theoretical flutter 
frequencies for stressed panels and thereby improve agreement between experimental 
and theoretical flutter boundaries in the transition region. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 9, 1965. 
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APPENDIX 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 

The International System of Units (SI) was adopted by the Eleventh General Confer- 
ence on Weights and Measures, Paris, October 1960, in Resolution No. 12 (ref. 16). Con- 
version factors for the units used herein are given in the following table: 



U.S. 

Conversion 


Physical quantity 

Customary 

factor 

SI Unit 


Unit 

(*) 


Length 

in. 

0.0254 

meters (m) 

j 

Pressure 

psf - lbf/ft 2 

47.88 

newtons per sq meter (n/tl 3 ) 

Stress 

psi = lbf/in 2 

6.895 X 10 3 

newtons per sq meter (N/m 2 ) 

Temperature 

(°F + 459.67) 

5/9 

degrees Kelvin (°K) 


* Multiply value given in U.S. Customary Unit by conversion factor to obtain 
equivalent value in SI Unit. 


Prefixes to indicate multiple of units are as follows: 


Prefix 

r 

Multiple 

giga (G) 

10 9 

kilo (k) 

10 3 

centi (c) 

10" 2 

milli (m) 

10 -3 
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TABLE I.- NATURAL FREQUENCIES MEASURED PRIOR TO EACH TEST FOR TEST PANELS MOUNTED IN PANEL HOLDER 


(a) a = 2.5. 
b 



fj, cps 

f 2 , cps 

f 3 , cps 

f 4 , cps 

q x 

Panel 1 

Clamped theory 

, 131 

154 

194 

258 

DO 


Simply supported theory 

/ 63 

90 

133 

194 

0 


Test 1 

j a 134 

158 

204 

273 

DO 


Test 2 

/ a 137 

160 

208 

276 

00 

Panel 2 

Clamped theory 

/ 147 

173 

219 

291 

OO 


Simply supported theory 

71 

101 

150 

219 

0 


Test 3 

132 

145 

198 

251 

62.5 


Test 4 

129 

140 

189 

246 

51.3 


Test 5 

a 149 

166 

227 

300 

OO 

Panel 3 

Clamped theory 

214 

251 

317 

422 

DO 


Simply supported theory 

103 

146 

218 

317 

0 


Test 6 

a 214 

242 

322 

419 

OO 


Test 7 

201 

225 

295 

393 

108.0 


Test 8 

203 

233 

307 

407 

129.0 


Test 9 

198 

231 

311 

415 

87.0 


Test 10 

197 

227 

309 

412 

85.2 

Panel 4 

Clamped theory 

214 

251 

317 

422 

OO 


Simply supported theory 

103 

146 

218 

317 

0 


Test 11 

204 

239 

316 

428 

143.0 


Test 12 

191 

227 

305 

412 

62.5 

Panel 5 

Clamped theory 

259 

304 

384 

512 

oo 


Simply supported theory 

125 

177 

264 

384 

0 


Test 13 

217 

250 

339 

462 

40.0 


o» l- 

= 2.9. 





fj, cps 

f 2 , cps 

f 3 , cps 

f 4 , cps 

q x 

Panel 1 

Clamped theory 

222 

249 

300 

376 

oo 


Simply supported theory 

105 

139 

214 

273 

0 


Test 1 

208 

231 

288 

370 

95.2 


Test 2 

183 

212 

278 

363 

36.4 


Test 3 

182 

211 

278 

361 

35.2 


Test 4 

179 

206 

277 

361 

33.8 


Test 5 

173 

203 

268 

353 

27.2 


Test 6 

173 

203 

268 

354 

27.2 

Panel 2 

Clamped theory 

145 

163 

196 

247 

oo 


Simply supported theory 

69 

91 

140 

179 

0 


Test 7 

a 143 

153 

196 

241 

570.0 


Test 8 

a 141 

102 

201 

— 

222.0 


Test 9 

a 145 

160 

199 

230 

OO 

Panel 3 

Clamped theory 

173 

194 

234 

294 

oo 


Simply supported theory 

82 

109 

167 

213 

0 


Test 10 

a 173 

190 

245 

313 

CO 


Test 11 

147 

162 

205 

279 

42.5 


Test 12 

154 

174 

224 

296 

57.2 


Test 13 

144 

161 

203 

273 

38.8 


Emitted in calculating average rotational restraint coefficient q x . 
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TABLE n.- PANEL FLUTTER TEST DATA, § = 2.5 

b 


[e = 10.5 X 10 6 psi (72.4 GN/m 2 ); a = 12.6 X 10' 6 °F' 1 ( 22.7 X 10" 6 


Test 

h 

T t 

q 

AT 

AP 

k \ 2 

aAT§| 

\hj 

|AP|/a\ 4 
E \hj 

2/3 

/ n \ 1/3 

( <1 \ a 

y3E J h 

* 

f, 

cps 

Flutter 
start 
or stop 

Panel 

condition 

(t) 

in. 

mm 

°F 

°K 

psf 

kN/m 2 

Op 

°K 

psi 

kN/m 2 

1 

0.063 

1.60 

310 

428 

2705 

129.5 

31 

17.2 

0.08 

0.55 

67 

35 

3.54 

72 

110 

Start 

F 

2 

.063 

1.60 

305 

425 

3910 

187.2 

17 

9.4 

.03 

.21 

37 

18 

4.01 

42 

125 

Start 

F 

3 

.071 

1.80 

305 

425 

2945 

141.0 

26 

14.4 

.26 

1.8 

44 

59 

3.25 

22 

135 

Start 

F 

i 





1610 

77.1 

78 

43.3 

.04 

.28 

133 

16 

2.66 

225 

100 

Stop 

B 






4090 

195.8 

118 

65.6 

.08 

.55 

190 

27 

3.61 

326 

300 

*Start 

B 

4 

.071 

1.80 ! 300 

422 

I960 

93.8 

16 

8.9 

-.04 

-.28 

27 

16 

2.84 

29 

120 

Start 

F 






1950 

93.4 

82 

45.6 

.03 

,21 

141 

13 

2.84 

236 

100 

1 Stop 

B 

5 

.071 

1.80 

300 

422 

1480 

70.9 

29 

16.1 

.16 

1.1 

50 

43 

2.59 

39 

115 ! Start 

F 






1480 

70.9 

75 

41.7 

.02 

.14 

129 

11 

2.59 

214 

105 

Stop 

B 






4650 

222.6 

126 

70.0 

.05 

.34 

217 

18 

3.79 

360 

290 

* Start 

B 






2505 

119.9 

152 

84.4 

.01 

.07 

261 

4 

3.08 

417 

310 

*Stop 

B 

6 

.103 

2.62 

305 

425 

3160 

151.3 

: 35 

19.4 

.01 

.07 

28 

3 

2.28 

1 42 

1 

170 

f Start 

F 

! 7 

.103 

2.62 

300 

422 

2445 

117.1 

89 i 49.4 

, .05 

.34 

i 74 

7 

2.09 

123 

— 

No flutter 

F-B 

1 8 

i .103 

2.62 1350 

450 

,2930 

! 140.3 

41 

22.8 

.01 

.07 

33 

3 

2.22 

1 49 

165 

Start 

! F 


! 




2920 

1 139.8 

54 

30.0 

1 .02 

.14 

43 

4 

2.22 

- 72 

155 

Stop 

B 

’ 9 

.103 

2.62 

355 

453 

2670 

1 127.8 

152 

84.4 

.01 

.07 

123 

3 

2.16 

j 197 


No flutter 

F-B 

10 

; .103 

2.62 

355 

'453 

4845 

232.0 

51 

28.3 

.05 

.34 

41 

7 

2.63 

; 58 

'170 

Start 

F 


! 




4830 

! 231.3 

1 

I 61 

1 

33.9 

.06 

.41 

49 

8 

2.63 

85 

160 

Stop 

B 

11 

.103 

2.62 

395 

475 

3640 

! 174.3 

42 

23.3 

| .03 

i .21 

33 

5 

2.39 

49 

160 

Start 

F 






2670 

i 127.8 

155 

1 

86.1 

| .03 

* .21 

125 

5 

1 2.15 

201 

130 

Stop 

B 

12 

.103 

2.62 

.485 

525 

1490 

: 71.3 

156 

86,7 

; .02 

i .14 

127 

1 4 

| 1.78 

1205 

1 

— 

No flutter 

F-B 

13 

.125 

3.18 

’490 

528 

4920 

235.6 

280 

155.6 

.05 

.34 

| 152 

4 

2.19 

244 

— 

■ No flutter 

F-B 


Small amplitude. 

t 

1 F denotes flat condition; B denotes buckled condition. 



TABLE HI.- PANEL FLUTTER TEST DATA, - = 2.9 

b 

(e = 10.5 X 10 6 psi (72.4 GN/m 2 ); a = 12.6 X 10~ 6 °F“ 1 (22.7 X 10" 6 °K" 1 J 

h T t <1 AT AP 9 r 71 2/3 i /o 

’est t — »atM iAPL/|\ (AS a , 

in. mm °F °K psf kN/m 2 °F °K psi kN/m 2 L ^ ' J v ' h 

1 0.081 2.06 350 450 4920 235.6 37 20.6 0.29 2.0 48 44 3.36 ! 

4865 232.9 121 67.2 .05 .34 158 14 3.35 2! 

2 .081 2.06 350 450 3905 187.0 43 23.9 .08 .55 56 19 3.11 I 

2475 118.5 132 73.3 .11 .76 171 22 2.67 21 

3 .081 2.06 375 464 2675 128.1 48 26.7 .05 .34 62 14 2,75 i 

2675 128.1 84 46.7 .04 .28 109 12 2.74 1' 

2890 138.4 101 56.1 .02 .14 131 6 2.82 2: 

2265 108.4 130 72.2 .05 .34 169 14 2.60 2‘ 

4 ,081 2.06 350 450 2755 131.9 123 68.3 .08 .55 160 19 2.77 2< 

3350 160.4 136 75.6 .27 1.9 176 42 2.96 2! 

5 .081 2.06 350 450 1970 94.3 118 65.6 .14 .97 150 27 2.50 21 

6 .081 2.06 370 461 2335 111.8 54 30.0 .09 .62 70 20 2.62 ll 

2330 111.6 60 33.3 .08 .55 77 19 2.62 l: 

2325 111.3 126 70.0 .10 .69 163 24 2.62 2i 

2500 119.7 130 72.2 .05 .34 169 14 2.68 2 


.053 1.35 315 430 2455 117.5 
.053 1.35 300 422 1680 80.4 


.053 1.35 305 425 3190 152.7 


2.2 .05 

9.4 .08 


17 9.4 .08 

18 10.0 -.01 

16 8.9 .36 


08 .55 

01 -.07 

36 2.5 


.063 1.60 310 428 3055 146.3 25 13.9 .11 

.063 1.60 300 422 3050 146.0 15 8.3 .10 

.063 1.60 310 428 1710 81.9 23 12.8 -.03 

1460 69.9 83 46.1 -.03 


83 46.1 -.03 -.21 

12750 I 131.7 1105 58.3 -.03 -.21 

3410 163.3 114 63.3 .14 .97 


.063 1,60 400 478 3675 176.0 


3655 175.0 130 72.2 .04 

4765 228.1 203 112.7 .00 


9.4 .15 1.0 

2.2 .04 .28 

2.7 .00 .00 


| i || | 3895 | 186 * 5 I 211 | 117 - 2 | * 02 | *14 

♦Small amplitude. 

F denotes flat condition; B denotes buckled condition. 
aQscillograph did not operate. 


/ q \ 1/3 a 
\p eJ h 

* 

f, 

cps 

Flutter 
start 
or stop 

Panel 

condition 

(t) 

3.36 

54 

150 

Start 

F 

3.35 

255 

95 

Stop 

B 

3.11 

80 

155 

Start 

F 

2.67 

280 

85 

Stop 

B 

2,75 

92 

145 

Start 

F 

2.74 

177 

130 

Stop 

B 

2.82 

210 

100 

Start 

B 

2.60 

273 

60 

Stop 

B 

2.77 

244 

90 

Start 

B 

2.96 

297 

70 

Stop 

B 

2.50 

280 

— 

No flutter 

F-B 

2.62 

102 

160 

Start 

F 

2.62 

131 

160 

Stop 

B 

2.62 

268 

80 

Start 

B 

2.68 

274 

90 

Stop 

B 

4.08 

52 

135 

Start 

F 

3.59 

80 

130 

Start 

F 

4.45 

9 

135 

Start 

F 

3.68 

64 

(a) 

Start 

F 

3.68 

32 

150 

Start 

F 

3.08 

72 

135 

Start 

F 

2.88 

290 

100 

Stop 

B 

3.56 

366 

110 

Start 

B 

3.82 

410 

110 

Stop 

B 

3.92 

35 

150 

Start 

F 

3.92 

450 

105 

Stop 

B 

4.28 

687 

270 

♦Start 

B 

4.00 

720 

270 

♦Stop 

B 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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